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The existence of structure on large (say, 100 Mpc) scales, and limits to 
anisotropies in the cosmic microwave background radiation (CMBR), have 
imperiled models of structure formation based solely upon the standard cold 
dark matter scenario. Novel scenarios, which may be compatible with large- 
scale structure and small CMBR anisotropies, invoke non-linear fluctuations 
in the density appearing after recombination, accomplished via the use of late- 
time phase transitions involving ultra-low mass scalar bosons. Here, we study 
the statistical mechanics of such phase transitions in several models involving 
naturally ultra-low mass pseudo-Nambu— Goldstone bosons (pNGB’s). These 
models can exhibit several interesting effects at high temperature, which we. 
argue are the most general possibilities for pNGB’s. 
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I. INTRODUCTION 


As the Universe is probed on larger scales, evidence for very large scale structures 
seems to be emerging. Voids, filaments, and walls on scales as large as 100 Mpc have been 
observed in various redshift surveys. 1 The existence of these structures, constraints on 
the anisotropies of the microwave background temperature, 2 and the existence of quasars 
at redshifts larger than 4, 3 makes it extremely difficult to understand the origin of large 
scale structure within the framework of the standard gravitational instability theory, i.e, 
cold dark matter with density perturbations coming from inflationary models. 4 

Many of the constraints on structure formation models can be obviated if a mecha- 
nism could be found which allows for density inhomogeneities to appear at a redshift z 
satisfying z^ c > z 1 with 8pj p ~ 1 (here z^c ~ 1000 is the redshift at recombination). 
In this case, the perturbations will not affect the CMBR directly (but may have signifi- 
cant indirect effects) yet structure on large scales will have ample time to grow so as to 
satisfy the constraint coming from quasar observations. 

Taking our cue from the fact that sources of density fluctuations may arise from the 
effects of phase transitions, it is interesting to ask whether phase transitions could occur 
at late times (i.e., after decoupling) in such a way as to generate large density fluctu- 
ations. This is not a new idea. Wasserman 5 showed that the existence of a first-order 
phase transition at late times could generate large fluctuations due to bubble collisions. 
Hill, Schramm, and Fry 6 proposed the idea of domain wall formation in late-time phase 
transitions in the context of ultra-low-mass pseudo-Nambu-Goldstone bosons that can 
readily occur in a wide class of models. 7,8 ’ 9 The implications of late time phase transitions 
with soft domain walls has been a subject of considerable activity in recent years. 10 

Press, Ryden, and Spergel considered the possibility of a slow-roll transition in a soft- 
boson model, 11 which not only drove structure formation, but also implied that the dark 


1 



matter is the residual oscillations of the field about the potential minimum. Schramm 
and Fuller 12 have also considered such scenarios within the context of Majoron models. 13 
These different approaches have a common theme: Ultra-low mass particles, typically 
spin-0 bosons, are a generic component of all such models. 

The most familiar spin-0 particles occuring in nature are the it mesons. The scale 
of the masses of the pions in comparison to the scale of nucleon mass is small, m 2 
m? N . This is well understood: the pions are pseudo-Nambu-Goldstone bosons (pNGB’s) 
associated with the dynamical breaking of fermionic chiral symmetries. In the limit of 
vanishing up and down quark masses, — ► 0, the pion masses go to zero, m 2 — ► 0, 

and the pions become exact Nambu-Goldstone bosons (NGB’s). Most of our intuition 
about pNGB’s derives from this established system, which is one of the most profound in 
elementary particle physics. We will exploit and develop the analogy with this system in 
greater detail in Section IV. The basic lesson is: Many continuous (perhaps approximate ) 
global symmetries may exist in nature that are spontaneously broken, and have associated 
NGB’s (or pNGB’s) with phenomenological implications. 

One such example is the familiar axion, 14 a hypothetical pNGB associated with the 
Peccei-Quinn (PQ) U(l) symmetry. The PQ-symmetry is broken by QCD instanton 
effects and the axion thus develops a small mass, m^a 0n ~ 0(m*f*/ f^don), where one 
conventionally assumes /«* ion ~ 10 12 GeV. Thus, the Compton wavelength of the axion 
is measured in centimeters, requiring the construction of macroscopic microwave cavities 
as detectors. Of course, since very low mass particles such as the axion are very difficult 
to detect, few theorists spend their time trying to invent new ones. However, the axion, 
which is a respectable, if not desirable, theoretical entity portends an important lesson: 
the physical world may contain many new phenomena in the far infrared which are not 
directly accessible, but nevertheless may play an important role in nature. 

A generalization of axions to a class of pNGB’s with masses of order m 2 ennion j f has 
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been analyzed in some detail. 7 Remarkably, if one associates mu union ~ Tra neu tnno ~ 0-01 
eV and / ~ M G ut ~ 10 15 GeV, one arrives at a cosmologically interesting scale A ~ 
f l m \ ~ Megaparsecs. This can lead naturally to a late-time phase transition. 

Embedded into theories that contain NGB’s or pNGB’s is our only rational guide- 
line in thinking about ultra-low mass particles: the principle of “naturalness.” In this 
regard, the mass scales of such particles must not be fine tuned, and must appear as 
a consequence of some plausible mechanism. 4 t Hooft first gave a concrete definition of 
the principle of naturalness: a parameter is “naturally ” small if when it is set to zero, 
the symmetry of the Lagrangian is increased. 15 In this case, the parameter will be multi- 
plicatively renormalized and will remain small to all orders of perturbation theory. While 
the cosmological applications are insensitive to whether or not a given model Lagrangian 
has been fine tuned, the form of any given low-energy effective Lagrangian, or its finite- 
temperature corrections, will be strongly influenced by the symmetries of the interactions 
of the full theory, and therefore we focus on natural models. 

More generally, there are two versions of the naturalness principle: (1) “Strong Nat- 
uralness,” in which the very low mass scales must emerge on the grounds of symmetry 
and dynamics without the input of any large hierarchy (for example, technicolor theo- 
ries respect this principle as a means of generating the hierarchy involving the W mass 
and the Planck mass, Mw/mpi ~ 10 -17 , although they have difficulty accomodating the 
observed large quark and lepton masses) (2) “Weak Naturalness,” in which one inputs 
a large hierarchy ab initio , which is then protected by a symmetry in the theory from 
being overturned by radiative corrections (for example, supersymmetry operates in this 
mode of protecting the hierarchy Mw/mpi , and “chiral” symmetries protect small ratios 
like TOneatimo /^electron ^ 10 )• 

The axion falls into the category of strong naturalness, since it would be an identi- 
cally massless particle by virtue of a symmetry principle if it were not for QCD effects 
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(instantons) which spoil the symmetry and are operant at energy scales of about 1 GeV 
(the QCD scale arises naturally from, e.g., any Grand Unified Theory upon specifying 
clqcd at the GUT scale). Other kinds of pNGB’s, having masses given by approximate 
expressions such as ~ Wermion//> where the decay constant / can be viewed as large, 
say 10 15 GeV to 10 19 GeV, are technically naturally low mass particles. 7 Here, the boson 
mass is protected by fermionic chiral symmetries such that if mf erm ion = 0> then 
vanishes to all orders of perturbation theory. 

Theories with naturally low-mass particles should be contrasted with theories where 
a small mass is unnatural. From a particle theorist’s point of view, the model of Press, 
et al., 11 suffers from being unnatural. The Lagrangian they considered assumes a mass 
term for a scalar field multiplet that is fine tuned to be of order (30 kpc) -1 , yet the 
field is assumed to have normal interactions with other particles. In any quantum field 
theoretic version of the model this would lead to an additive quadratic divergence in the 
mass term. Thus, to maintain the small mass term one must fine tune the theory in each 
order of perturbation theory. 

As a general laboratory for the statistical mechanical phenomenology of pNGB’s, we 
will focus on the models developed by Hill and Ross. 7 These models have a light pNGB 
^ which couples to fermions. The effect of these fermions is to induce a potential for <f>, 
which can lead to a phase transition. These models are very simple, but we believe that 
they are sufficiently general to imitate any kind of pNGB dynamics. For example, the 
Zff models for N > 2 lead to a phase transition analogous to the axion case. We remark 
that in this analysis we will not include the potential effects of anomalies, aside from 
briefly indicating in Section IV how they arise in pNGB physics. 

The interesting setting for these theories in a cosmological context is one in which 
the fermions are the light neutrinos and <f> is a NGB associated with symmetries of the 
neutrino masses. We will discuss this below in the context of pNGB’s associated purely 
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with hypothetical Dirac mass terms, as well as those associated with Majorana mass 
terms . 13 ’ 18 The critical temperature of the transition in some cases is either naturally of 
the order of, or determined by, the small masses of the neutrinos, and as such would 
be rather small compared to the usual scale of critical temperatures in particle physics 
models. 

The purpose of this paper is to understand the statistical mechanics of pNGB phase 
transitions. As a “straw man” the first and simplest case we present in Sec. II is an 
unrealistic one. The model presented in Sec. Ill is more realistic but technically unnatu- 
ral; it does however illustrate several features that will be present in more sophisticated 
models. It is analogous to a Coleman- Weinberg 17 effective potential with thermal or 
finite density corrections coming from relic neutrinos, and it has some striking features 
in common with slow-roll inflationary schemes. Sec. IV contains a discussion of the mo- 
tivation for neutrino pNGB models. In Sec. V, we review the model 8 where we give a 
standard computation of the effective potential in the tadpole formalism (which we use 
throughtout ). 18 We then study the finite temperature effects for this model. The usual 
trick of using the high temperature expansion in order to determine Tq is unreliable and 
more delicate methods must be used. We then generalize the Zi model to Zff models, 
which softens the fermion loop effects in the UV. We find that these models do not un- 
dergo a conventional phase transition. The potential “turns on” at low temperatures in 
analogy to the axion case. We summarize our results in the final section. 


II. SELF-INTERACTING SCALAR MODEL 


Let us review a phase transition associated with a simple model consisting of a single 
real scalar field with a classical potential of the form 
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( 2 . 1 ) 


W) = + j-W*. 

The classical potential Eq. (2.1) has minima at <f> — ±<r, where <r = ^Jm^/X 0 . The mass 
of the scalar field is related to the curvature of the potential at the minimum: 


m 


, = ^Vq(^) 
*- ip 



= 2 ml = 2A 0 <r 1 . 


( 2 . 2 ) 


This is a very simple model that illustrates the phenomena of phase transitions and do- 
main wall production. The calculation of the critical temperature of the phase transition 
in this model is well known and completely straightforward. We review it here to estab- 
lish some notation and definitions that will be of use in the more complicated models 
discussed below. 

Questions of symmetry breaking, symmetry restoration, finite-temperature effects, 
etc., are best studied by considering the “effective potential.” This will account for the 
quantum effects of virtual particle emission and absorption, as well as the effect of emis- 
sion and absorption of particles from the thermal background. Methods of calculating 
the effective potential are well developed. In one prescription the evaluation of the po- 
tential involves shifting the field by an arbitrary amount (say <f> — » <}> + $), and evaluating 
the “tadpole” diagram of Fig. la in the shifted theory. In this formalism the effective 
potential to one loop is 


v{4>) = v 0 (4)- [d$iw\ , 


(2.3) 


where I'M is simply a factor of i times the tadpole diagram of Fig. la in the shifted 
theory. In the shifted theory the potential is 


V 0 (<l>) = ~m 2 0 (<f)- 4>f + $)*, (2.4) 

which results in a coupling constant for the cubic term of A <f> and a mass-squared of 
— ttIq + 3A0 2 . Evaluating the tadpole of Fig. la, I'M is simply 
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Integrating with respect to <f> and rotating to Euclidean momentum k y the one-loop 
correction to the potential is 


W) = i/^jK* 1 -m5 + 3A^ ) . (2.6) 

This expression of course is divergent and it is necessary to cut off the integral at A: 2 = A 2 , 
with result 


v{4>) = W) + W) 


— — Tn^<f> 2 + — A<£ 4 + a<f> 2 + b<f>* + 


( —m\ + 3A<£ 2 ) 2 ( —ml + 3A<ft 


64tt 2 


in(: 


A 1 


■(2-7) 


where a and b include terms proportional to A. These constants will be determined 
by renormalization conditions, e.g., by the definition of a renormalized mass m and a 
renormalized coupling constant A. After renormalization, the zero-temperature, one-loop 
potential may be expressed as 


Vi(<« = 


i 

647T 2 


In 



( 2 . 8 ) 


where ft is an arbitrary mass scale which can be related to the renormalized coupling 
constants, and M 2 (<f>) = —ml + 3A<£ 2 is the mass as a function of <f> . 

The finite-temperature 19 corrections to the potential arise from the interaction of the 
(f> field with the ambient background. To calculate the effect of the background, one 
computes the quantum corrections to the tree-level potential of Eq. (2.1), taking into 
account the fact that the background influences the <f> propagator. That the background 
should have an effect at the one-loop level is easy to see, since evaluation of the effective 
potential in the one-loop approximation involves evaluation of the tadpole diagram of 
Fig. la, which in turn involves the <f> propagator. The <j> propagator is influenced by the 
distribution of real particles in the background. 



If the phase space density of the <j >' s is denoted by f<f,(k), the <f> propagator (in the 
real-time formalism) becomes 

D T (k) = i(k 2 - M 2 )- 1 + 2irU(k)6{k 3 - M 2 ). (2.9) 

Consider the part of the momentum integration of the tadpole diagram for emission and 
absorption of a particle on shell ( k 2 = M 2 ). There is no way to differentiate between 
the possibility that the absorbed particle is the virtual particle emitted, or the absorbed 
particle comes from the background. The second term in Eq. (2.9) accounts for the latter 
possibility. 

Clearly the effect of the background particles depends upon their phase-space density. 
If the <^’s are in thermal equilibrium, they will be distributed in phase space according 
to the Bose-Einstein distribution: = [exp(J5/T) — l] -1 , where E = y|kp + M 2 . 

For the moment, we will make the assumption that the phase-space distribution of the 
^’s are described by the equilibrium expression. 

In the one-loop approximation, the potential is a sum of the tree-level potential, Vo(^) 
given by given by Eq. (2.1), a zero-temperature one-loop correction, of Eq. (2.8), 

and the temperature-dependent one-loop potential, A Vy(<^>): 20 

V(<£) = Vo (0) + W) + A V t (4>). (2.10) 

The temperature- dependent part of the propagator adds to I^ 1 ) a term 

r£> = -\J ^;6 Xj, 2ir U(k)6[k' - (- ml + 3A?)]. (2.11) 

Following through the integration with respect to rotation to Euclidean momentum 
A, and integration over d 4 fc, one obtains the (finite) result 

AW) = ^ jT <***’ln [l - exp (-^Ttf*?)] , (2.12) 
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where again, M 2 ((f>) = — m\ + 3A^ 2 . 

To demonstrate that there is a phase transition and to calculate the critical tempera- 
ture is straightforward. At zero-temperature the minima of the potential are <j> = ±<r , and 
the curvature at <f> = 0 is negative (i.e., <j> — 0 is a local maximum). At high temperature 
AVy(^) can be expanded in <f>, with the leading-order ^-dependent term proportional to 
+T 2 <f> 2 . Clearly at high temperature the curvature of the potential at (f> = 0 is positive, 
and indeed <j> = 0 is the true minimum of the theory at high temperature. We will denote 
the temperature at which the curvature of the high-temperature minimum vanishes as 
the critical temperature, T c . In the above theory, d 2 V/d<f> 2 evaluated at <)> = 0 changes 
sign at a temperature Tc = 2 <7. 

This model illustrates the standard scenario for making walls. At temperatures above 
the critial temperature, the value of the field is pinned at the high-temperature minimum, 
<f> = 0. This is because at high temperatures <f> = 0 is the global minimum of the potential, 
and furthermore, the mass of the field at high temperature is large (of order AT). It is 
the large mass that pins <f> to the high-temperature minimum. Now once the temperature 
drops below the critical temparature, the <f> field will evolve classically to either of two 
possible minima. Regions of the Universe in different minima will be separated by a 
domain wall. This scenario depends upon the fact that as the phase transition starts, <j> 
is localized at a low- temperature maximum, which is also a high-temperature minimum. 

Note that Tc/m^, = \J2/ A. Thus, it appears that by making A sufficiently small, one 
mi gh t have a late-time transition generating soft-walls. Let us explore this scenario. 

Let us assume a generous range for Tc, say To < Tc < T rec , where To is the present 
temperature, T 0 = 2.7 K ~ 2.4 x 10 -4 eV, and T rec is the temperature at recombination, 
Tec ~ 0.3 eV. Let us also assume that the boson has an ultra-low mass, < 10 -24 eV. 
The combination Tc ^ T and m $ < 10 _24 eV leads to the constraint A < 10 -4 °. If this 
constraint is satisfied, the model as presented will lead to a late-time phase transition, 
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and contains soft domain walls with thicknesses of order m^ 1 ~ parsecs. 

However, there are two serious problems with the model. The first problem is that 
it is unreasonable to assume that a fundamental constant, such as A, has a value of 
10 -40 without some deeper underlying motivation. Such a value is unnatural in the 
technical sense 15 and is arbitrary. A second difficulty is that with such a small value 
of A it is unlikely that the 0’s were ever in equilibrium and the assumption that they 
are present in a thermal phase-space density cannot be justified. Unless the potential of 
Eq. (2.1) is augmented by some additional interaction terms, the only processes leading 
to thermalization of the <f> y s are <f> self-interactions. The cross section for this processes in 
the relativistic limit is <r m t = A 2 j a ~ A 2 /T 2 , where the last approximation assumed that 
the average energy of the (j> is characterized by a temperature T. If <f> is in equilibrium 
and relativistic, n $ ~ T 3 , so the interaction rate of the <£’s is T xn t ~ ~ A 2 T. In 

the radiation- dominated era, the expansion rate is H ~ T 2 /mph so T- m t/H ~ A 2 mpijT. 
If A = 10~ 40 , then T^t/H < 1 for T > 10 ^mpi < T 0 . A similar conclusion follows for 
the expansion rate appropriate to a matter-dominated era. 

Clearly the assumption that <£’ s exist as a thermal background cannot be justi- 
fied. Of course, one might imagine that the background is not established through 
self-interactions, but rather is the result of some non-standard (but reasonable) process 
such as primordial black-hole evaporation, quantum effects during inflation, or other such 
processes. 

We will now turn to our attention to developing models where the phase transition 
is driven not through <j> self interactions with a background, but rather by <j> interactions 
with a background of some other field “0, typically a fermion. The virtue of this compli- 
cation is that it is possible to have ^ interactions weak enough to provide a late-time, 
soft- wall transition, but the if? can have additional interactions that can establish the 
background by thermal interactions. 
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III. SCALAR FIELDS WITH YUKAWA INTERACTIONS 


To the classical potential of Eq. (2.1) we add a Yukawa coupling of ^ to a fermion 
field ip: 31 

Vo(<f>,ip) = + ^<t>* ~ (3.1) 


where the parameters mo, Ao, and h are the unrenormalized mass and coupling constants. 22 

Before turning to the temperature-dependent effects, consider the zero-temperature 
radiative corrections. The one-loop corrections involves calculation of the tadpole dia- 
gram of Fig. lb in addition to the scalar tadpole of Fig. la. In the following we will 
assume that the fermion loops dominate, which will be true if if h 2 Ao, and ignore the 
boson tadpole. The effect of the <p-xp interaction on the effective potential is evaluated 
by calculating the tadpole as discussed in the previous section. Upon shifting the field 
<f> — ► <f> + the mass of the %p is M ^ = h<p> and the <p-Tp-Tf> vertex that appears in the 

tadpole is proportional to h. Thus I'M is obtained from computing the one-loop tadpole 
diagram of Fig. lb: 


r (U = i f (3.2) 

J (2tt) 4 £ - h<f> 

Following a procedure similar to the one outlined in the previous section, one obtains 
terms in the one-loop effective potential that are infinite (proportional to a cut off A) 
and terms that are finite. The infinite terms are dealt with by some renormalization 
prescription, and the renormalized, one-loop effective potential is simply 


v,{4>) = 


1 (M 4ln ^> 


167T 2 


(3-3) 


where fi again is an arbitrary mass scale related to the values of the coupling constants. 

Now the effect of a background of real on the effective potential is calculated along 
the line as the previous section. Again the tadpole of Fig. lb is calculated replacing the 
fermion propagator by its finite-temperature expression: 
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S T (k) = i()S - M t y' - 2*U(k)(l + - Ml), 


(3.4) 


where here is the phase space density for xj). Again let us assume the phase space 
density for xp is that of a thermal distribution (i.e., a Fermi-Dirac distribution with 
temperature T). This adds to the one-loop effective potential a temperature-dependent 
term 

AVr(<£) = -4 ^2 J Q d x x2 ^ n i 1 + ex P ir^j x2 + M H T2 )\ • ( 3 * 5 ) 

In comparison to Eq. (2.12) several differences are obvious. The overall sign is opposite 
because there is an overall sign difference between fermion and boson loops. The sign 
difference in the argument of the logarithm arises from the sign difference in the Fermi- 
Dirac verses Bose-Einstein distribution functions. Finally the overall factor of 4 owes to 
the trace over 7-matrices involved in the fermion loop. 

Now let’s consider the phase transition. Let us assume that the zero- temperature 
one-loop potential has negligible effect and the curvature at <f> = 0 remains negative. 
Expanding Eq. (3.5) for large T, A Vj(<£) ~ -{-h 2 <f) 2 T 2 /3. Clearly at high temperature 
the curvature of the potential at <f> = 0 will be positive, and again it will be the global 
minimum of the potential. We again define the critical temperature for the phase tran- 
sition to be the temperature where d 2 V/d<f> 2 evaluated at <j> = 0 vanishes. This results in 
a critical temperature of T c /m + = 3. This expression is very similar to the critial 

temperature in the model of the previous section with the replacement h +-+ \/X 

Clearly by making h sufficiently small it is possible to have Tc for a late- 

time, soft- wall phase transition. However the present model is superior to the previous 
model in one important regard: Although the xp field driving the transition must be very 
weakly coupled to (p y it may have stronger couplings to other fields. These other (yet 
unspecified) couplings can be sufficiently strong to establish xp in thermal equilibrium. 
Therefore, although <p may be completly decoupled from the thermal bath, its interactions 
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with the thermal background of ip's can restore the symmetry at high temperature and 
lead to a phase transition. 

The model still suffers from the ugly feature of very small, unnatural dimensionless 
coupling constants. In the next section we will describe the physical motivation for the 
origin of such small dimensionless coupling constants. The model will have some of the 
features of the model of this section. Before proceeding, let us restate the parameters of 
the present model. The model has a scalar field <f> with mass m $ < 10~ 24 eV. The vacuum 
expectation value of the scalar field is a. The Yukawa coupling of 0 to a fermion field 
ip results in a mass = her. If h is much larger than the scalar quartic self coupling, 
the phase transition temperature will be Tc ~ m$fh. If we want Tc to be larger than 
To , then h < 10 20 . However we are free to choose ar to be as large as desirable, and 
to have the ^ field coupled to other particles with sufficient strength to establish it in 
equilibrium. 


IV. PSEUDO-NAMBU-GOLDSTONE BOSONS 
A. Chiral Lagrangians 

The model discussed in Sec. Ill most simply demonstrates the basic idea of a late-time 
phase transition, but it suffers from the lack of symmetries that can naturally give a soft 
boson mass scale for <f> without fine tuning. Let us now consider models in which these 
constraints are implemented. 23 

Consider first the low-energy effective Lagrangian which contains a neutrino field v. 
C = + VLi<f)v L + Vri^vr + ( mVit/Re + h.c.) (4.1) 

4U 

where ui (i /r) is the left-handed (right-handed) projection: */£ = (!— 7 s)i //2 (i/r = 
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( 1 + 75 )*'/2)- The factor of me 1 *!* can be viewed as arising from the vacuum expectation 
value (VEV) of some Cf(l) complex scalar field $ that is coupled as gvLVR$ + h.c. In a 
17(1) invariant potential V($) we assume that $ develops a VEV of ($) = /e*^/V 2, 
and m = gf/V% (the factor y/2 assures that (f> has a properly normalized kinetic term of 
( 50) 3 /2 coming from the kinetic term of $, |#$| 3 ). 

Eq. (4.1) is a “chiral Lagrangian,” possessing the continuous chiral 17(1) symmetry: 

vl — * e lOL UL\ vr e~ ta i/R\ <j> —* <f> + 2a/. (4.2) 


We emphasize that the symmetry is not broken, and is properly said to be “nonlinearly 
realized” (this is often a confusing point: spontaneously broken symmetries are in fact 
equivalent to nonlinearly realized symmetries and are not really broken symmetries). We 
remark that chiral Lagrangians have several important and well-known properties: (1) 
as stated above, they can be embedded into a fully renormalizeable theory in which, e.g., 
a t 7 (l) complex field develops a vacuum expectation value, ($) = f/y/ 2 , and <f> is then 
the residual Nambu-Goldstone boson; (2) C can itself be viewed as renormalizeable for a 
small cut off A <C / up to suppressed counterterms of order A//; (3) <j> will be identically 
massless unless terms are introduced which explicitly break the chiral symmetry; (4) (f> 
satisfies “Adler decoupling,” i.e., we may replace v everywhere by v': 


u L e»'”-, 


»' R = i *e-W; 


(4.3) 


and our Lagrangian becomes: 

C = + V> L i$ v' L + v' R + m(v' L v' R + h.c.) + id'V (4.4) 

and we thus see that <j> disappears in the mass term but couples derivatively to the 
neutrino as d^cf) 1 / 757 ^ 1 /. Therefore, for small <j> momentum <f> emission or absorbtion 
amplitudes will tend to zero as q^ — ► 0 (“Adler zero”). An implication of this is that 
<j) will not mediate a long-range 1 /r 2 force as a consequence of this decoupling theorem 
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(when the symmetry is broken by a nonchirally invariant mass term as discussed below, 
then the Adler decoupling can be violated and <f> can mediate a long range force, though 
this requires CP- violation 7 ). 

We remark that at this stage on yet another branch of NGB or pNGB physics. If 
we integrate the last term of Eq. (4.4) by parts we obtain — ( 2 /) _1 <£ ^ 1 / 757 ^ 1 ^. Now if 
we include the effects of generic gauge fields that may be coupled to v (for example, the 
vl couples to the electroweak gauge fields) then the divergence of the axial current will 
contain an axial anomaly ~ cF F + ..., and therefore we find that <f> couples 

to the gauge fields through this anomalous term: — (2/)“ 1 0 cFF . This is an explicit 
symmetry breaking effect coming from quantum loops and it generally leads to important 
consequences. For example, the decay 7r° — + 27 involves this term; the gluon field enters 
the divergence of the PQ current ultimately giving a mixing of the axion to the 7r°, tj and 
77 ' from which the axion mass derives. In principle we should include potential anomaly 
effects in our effective potential analysis, however we will not do so for a reason: the phase 
transitions we consider here occur at very low temperature (or finite density) and arise 
from other explicit symmetry breaking effects. It is hard to see how anything but the 
electromagnetic anomaly could play a role at these low energy scales. It is conceiveable 
that an electromagnetic effect, e.g., in a plasma, might trigger a late-time transition 
through the anomaly, but we will not consider this possibility in the present paper. 

Let us now consider the explicit breaking of the symmetry by effects other than 
anomalies. By this we mean the addition of new terms to Eq. (4.1) which explicitly 
violate the nonlinearly realized symmetry of Eq. (4.2). For example, to the Lagrangian 
we may add a small mass term for <f> of unspecified origin. Usually this comes from some 
deeper symmetry breaking in the theory which breaks the continuous U( 1) down to a 
discrete subgroup Zn. For example, let us break U(l) to its trivial center by adding a 
“soft-breaking” term, which is a cosine potential for <j > . This implies that ($> — ► <p + 2mr f 
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remains an invariance. So we now have: 


£' = £ + m* cos + 0). 


(4.5) 


By expanding the cos <£// term about a local minimum we infer the mass of the (j) boson: 

= to 4 // 2 ( 4 - 6 ) 

and there are also further interaction terms such as a \<f > 4 term where 


A = m 4 /12/ 4 . (4.7) 

The physical values of A and are proportional to the ratio m//, and can be almost 
arbitrarily small, while remaining stable under quantum radiative or thermal corrections. 
This is “natural” in the sense of ‘t Hooft, 15 and is due to the fact that m = 0 is a symmetry 
limit of the full theory [in which we recover Eq. (4.2)]. 

For example, with / ~ 10 16 GeV and m ~ 10 -3 eV we have m# ~ 10 -3O eV, or 
a Compton wavelength, h/m^c ~ 10 Mpc. The incoherent particle interaction rates 
will be negligible since A ~ 10 -109 ! The Adler decoupling theorem still holds with soft 
breaking since it follows from redefinition of fermion fields. In general, reaction rates 
involving <f> coupled incoherently to matter will be suppressed, since the cross sections 
are necessarily proportional to a power of 1 // 3 . Thus it is difficult, if not impossible, 
to excite <f> in the laboratory, just as the detection of invisible axions is difficult. Since 
reaction rates that maintain thermal equilibrium of <f> are of order T 3 // 3 , in a radiation 
dominated Robertson— Walker phase we see that the condition that <f> be in equilibrium 
is T 3 // 3 > T 3 / mpi or T > f 2 /mpi. Hence, a pure pNGB like (/> decouples very early in 
the evolution of the Universe. 

What kind of deeper structure can give rise to a mass term for <f>l In the case of QCD 
the proton and neutron are analogues of the v field, and the pion is the analogue of <f>. The 
deeper structure that breaks the chiral symmetry is the presence of light quark masses, 
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which are not chirally invariant. This leads to the non-zero pion mass term. However, it 
is unlikely that the only manifestation of a deeper symmetry breaking term is merely a 
mass term for <f) . Indeed, in the case of the nucleon-pion system, the finite quark masses 
also lead to a small chiral symmetry breaking term in the proton and neutron masses 
(known as the {r-term). We can make a strict analogy to this situation in the present case 
by adding an additional neutrino mass term to the Lagrangian, which explicitly breaks 
the chiral symmetry , in analogy to the QCD cr-term. The low-energy Lagrangian then 
becomes: 

C - ^d^cf>d^(f> + V Lift v L + VRi<f) ur 
z 

+(rrwLVRe lli> ^ + €vlVr + h.c.) + m* cos (<f>/f + 6 ), (4.8) 

where the term involving e explicitly breaks the symmetry of Eq. (4.2). Now, if m — * 0 
we must also set e — ♦ 0 to recover the symmetry limit of Eq. (4.2). However, a nonzero to 
will always be induced by the presence of a nonzero e and m. For instance, the diagram 
of Fig. 2a with a cutoff A < / gives an induced term in the Lagrangian 

Anduced ~ ~^~2 COS (^//)* ( 4>9 ) 

In the present case we see that the induced scalar mass will be of order: 

~ me(A 3 // 3 ) ~ me. (4-10) 

We can view this as the origin of the scale of m 3 ~ >Jme A. The mass can be naturally 
small in the technical sense since we can tune the symmetry breaking parameter e to be 
arbitrarily tiny for large m so that the observed neutrino mass is, e.g., m„ ~ mo ~ 1 eV, 
while m<j> ~ (100 Mpc)"* 1 with 6 ~ 10~ 60 eV, and the symmetry will guarantee that we 
don’t have to worry about radiative corrections changing this result. This is arbitrary, 
however, and this is not the ultra-low mass case we seek for application to a late-time 
phase transition. 
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In the Lagrangian of Eq. (4.9) we observed the appearance of a (“large”) quadratically 
divergent contribution to the induced mass of <j>. Can we somehow reduce the degree of 
divergence of this induced term? The answer is yes: residual symmetries can readily 
control this. 

Consider the following Lagrangian containing N Dirac neutrino species and invariant 
under a Zn discrete symmetry: 

C = -d^<j>d u (f> + T. Vjifii'j + f m + VjL^jR + h.c. (4.11) 

2 j= 0 j=0 

The continuous U( 1 ) chiral symmetry is broken down to a residual Zy discrete symmetry: 

Vj — ► u i+ 1; i/N-i -* u 0 \ </> — ♦ 4> + 2ir f/N. ( 4 -12) 

If one now computes the induced <p mass term, one obtains the (^-dependent term 

N-l 

(*•“) 

j=0 1U7r 

where MJ = m 2 + e 2 + 2mecos(<f>/ f + 2jir/N ). Notice that the potential retains the 
discrete symmetry <j> — > <j> + 2j t Kf/N. Now, it is readily seen that is a constant 

independent of <f> for N > 2. Therefore the A dependence in Eq. (4.13) is illusory; the 
^-dependent part is A independent, and for N > 2 we may write 

N-l JufA 

v {4>) = - J2 log ( M i) + const. (4.14) 

Hence, in Zs models we can view the symmetry breaking as soft and the potential of <f> 
is calculable. 

These models can be further generalized. In Ref. 7 the effects of CP violation are 
also included to contruct models in which the Adler decoupling theorem is violated and 
the pseudo-Nambu-Goldstone bosons develop CP-violating Yukawa couplings. This is 
analogous to including a 0-term into QCD (without an axion to kill it!). The net effect is 
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the possibility of weak, sub-gravitational strength long-range forces in a natural model. 
This is a further complication of the model which we will not include at present. 


B. Majoron Models 


In discussing neutrino masses the most sensible framework is that of the “seesaw” 
mechanism. 18 Here one attempts to explain the comparatively tiny values of neutrino 
masses relative to their charged lepton counterparts by invoking ultra-large Majorana 
masses for the right-handed neutrinos. In short, one needs only to assume that (i) all 
neutrinos have Dirac mass terms (perhaps of the same order as their charged lepton 
partners within a given generation), and (ii) all right-handed neutrinos have a Majo- 
rana mass term. With both terms present we have a conventional Gell-Mann-Ramond- 
Slansky-Yanagida 16 see-saw mechanism. The predicted light mass scale of neutrinos will 
be acceptably small, of order ^ 2 cpt on/^- The right-handed neutrinos are favored for a 
large Majorana mass term because they carry no known gauge symmetries. Nonetheless, 
one can also invoke small left-handed Majorana masses, as in the Gelmini-Roncadelli 
model. 24 These carry electroweak isospin of I = 1 and must be very small, since 7 = 1 
effects are suppressed in the Standard Model. 

Here we will give only a brief toy model discussion, leaving a more detailed catalogue 
of schemes to another place. 25 We consider a single Dirac neutrino field, with left- and 
right-handed components vl, Ur. We now assume the existence of Dirac and both left- 
and right-handed Majorana mass terms: 


-( UL \ ( e 771 
2 \ u r / \ m Me'Mf 



( 4 . 15 ) 


Here superscript- C denotes charge conjugation. We assume e <C m <C M. The phase 
exp(i^//) is the CMP Majoron, 13 the NGB associated with spontaneous breaking of the 
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rh-neutrino number symmetry. We take e to be an explicit U{1)l lh-neutrino 
number breaking effect. The Dirac masses communicate the explicit breaking in the rh 
sector into the lh sector, and thus the majoron becomes a pNGB. 

If we set either m = 0 or e = 0 then the phase exp (i<f>/ f) can be eliminated from 
the mass matrix by a redefinition of the neutrino fields, and we are left only with the 
derivative coupling and <j> remains massless. This is the usual assumption for the majoron. 
However, we see that the diagram of Fig. 2b implies an induced mass term for <f> given 
by 

em 2 M log(A 2 /m 2 ) /. 

Anduced COS (<^//)» ( 4 * 16 ) 

which is highly suppressed owing to the combination of small e and the chiral suppres- 
sion involving m 2 . Thus, we expect that majorons will behave in a mode which is no 
more divergent than the Z 2 case described above, and in more general schemes will be 
suppressed as in the Z ^> 2 case. 

In fact, the full behavior of broken majoron models may be very rich. Bjorken (pri- 
vate communication) has suggested considering the full 3-generation standard model to 
contain a spontaneously broken SU(3)r (which will be in a sextet mode, corresponding 
to v' r Vr) which will produce NGB’s. The SU{Z)r is then explicitly broken by the Dirac 
mass terms, and no lh majorana masses are included. The result is a hierarchy of NGB’s, 
some remaining massless while others acquire a spectrum of induced mass terms. More 
general schemes such as this will be considered elsewhere. 25 

We turn now to the thermal corrections to the chiral Lagrangians. 
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V. THERMAL PROPERTIES OF Z N MODELS 


In this section we will analyze the thermal properties of the Zw symmetric chiral 
Lagrangian models discussed in the previous section. For simplicity we will first consider 
the Z 2 model, and then generalize our results to the Zn case. 

A. Z 2 -symmetric models 

The effective low-energy theory (here, low-energy refers to scales much smaller than 

/) consists of two fermions (presumably neutrinos) j = 0,1, coupled to the scalar 

field <f> by Yukawa couplings of the form: 

1 

-£yuk = '52i’j{m + e[cos(4>/f + jir) + ij 5 sin(<f>/ f + jir)}} (5.1) 

j=o 

where we have used ^£,V'fie‘“ + h.c. = ify cos(a)+t'075'0 sin(a). Here, / can be thought of 
as the scale at which the continous symmetry, of which <j> is the Nambu-Goldstone boson, 
is spontaneously broken. The origin of and motivation for considering such theories was 
discussed in Sec. IV. 

We can rewrite Eq. (5.1) by performing a chiral rotation to eliminate the 75 term, 
with result 

— £yuk = — M + (<j>)Tj} 0 Tpo — 

M±(<f>) = m 2 + e 2 ± 2mecos(<^//). (5-2) 

Note that at this point there is no potential for <j>. However the effect of the coupling 
will generate a non-trivial potential for <f> through radiative corrections, rendering it a 
pseudo-Nambu-Goldstone boson. We first calculate the zero-temperature potential, then 
consider the finite-temperature potential. We will employ the same methods developed 
in Sec. III. 
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We will only consider the fermion contributions.^® The tadpole method described in 
Sec. II can be adapted to the present case. Rather than defining tadpoles in <f>, i.e., (<f>), 
as in the self-interacting scalar model, here we must preserve the full symmetry of the 
theory and define the tadpole to be expectation value of the mass terms. The one-loop 
potential now receives contribution from two tadpoles with fermions as in Fig. lb, from 
the mass terms (M+(^)^ 0 V» o) and These contribute to the unrenormalized 

one-loop potential a result given by 


W)= E 

i=+.- L 


A 4 


32tt 2 


87T 2 


- i) • < 5 - 3 ) 

Here A is an ultraviolet cutoff for the theory; presumably it cannot be larger than /, since 
at this scale the effective theory with the Nambu-Goldstone boson <f> must be supplanted 
by the full theory. 

We must now renormalize the potential. In so doing we will introduce an arbitrary 
energy scale fi\ of course no physical effects will depend on fi. To proceed, first introduce 
the scale /z into the potential: 


VM) = £ 

j= + t~ L 


A 4 




32ir 2 8 tt 2 




(5.4) 


Before proceeding we make note of the following identities: 


£ MjM = 2 (m’t! 1 ) 

J=+»- 

£ = 2[(m 1 + e 1 ) 1 + (2me)‘cos 1 W/)]. (5.5) 

j=+ 1- 

Because of the residual Zj symmetry is independent of <f>. From Eq. (5.4) 

we see that we must add counterterms Vct(<A) = Vo + Vi cos 2 (<^//), where Vo and Vi are 
(j) and n independent, to the original Lagrangian to cancel the cutoff-dependent terms. 
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Here V 0 and Vi are given by 

Vo = + + *)'*£ 

3 i » — 

V, = M‘( / *) + i i;2(2m e ) I ln^ ) (5.6) 

with Vj,(/*)i and M 4 (fi) finite. The final potential is then given by 
V{4>) = Vct(<^) + ^i(^) 

= V 0 (n) + cos 2 (d>//) 

-Tib [«*> ('" ^ - 1 ) + ^ - i)] • (5 - 7) 

Now the ^-independent parts of Vq(p) and M A (p) can be fixed by renormalization 
conditions. In particular, let us choose the renormalization conditions 


v"{4>%= o = m oJ 


F(tt/2) = 0. 


At this stage the sign of m* is not fixed. The final potential becomes 

T r,is 1 , 3 , 2x2 b 1^ 

v(« = 32iij("* + + m -) (k— V 2) 

+ -\f™l + ^i( m l - m -> ( m + ln ^r - m - ln 7r)] cos, W^ 


-lb H ( ln ^ - 5) + M -« ( h ^ - 5)] • (5 ' 9) 

where we have made yet another definition, m± = (mi e). Althought it is not apparent, 
Eq. (5.9) is p- independent, as can be seen by showing V(<j > ; p) — p') = 0. We leave 

the exercise in algebra to the reader. 

Obviously V(<f>) is periodic with period 7r and that its extrema are at <f> = 0, 7r/2 

(mod tt). The location of the minima depend on the sign of m\. We show the potential 
in Fig. 3 for negative m^. 

Let us now turn to the finite temperature corrections to the effective potential for <f>. 
Given the (^-dependent masses and we can then use the finite-temperature 

formalism discussed in Sec. Ill to compute the corrections [cf., Eq. (3.5)]: 
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AVr(<A) = - 4 ;tt X) / dx x3ln l 1 + ex P (-\/* 3 + Af i(^)/ 72 )] • 

27T 4 . . Jo L 

J — T 1 


(5.10) 




The signal of a second order phase transition is the flattening of the potential at the 

high-temperature minima, i.e., V"(<f> = 0)| r _ Tc = 0. Here, V(<f > ) = + A Vr{<f>). The 

temperature- dependent mass squared at <f> = 0, m 2 (T), is given by: 

4 meT 2 ^ /■« , , (x 2 + m^/T 2 )- 1 / 2 

+ exp(z 2 + m^/T 2 ) 1 / 2 

Since 4me = — ml > 0, it is easy to show that the temperature-dependent term 

is always positive . Thus, if <f> = 0 is a minimum at zero temperature, it will remain so at 
any finite temperature. This implies that the T = 0 maximum at 7r/2 (when ml > 0) 
remains one at finite T. Thus we do not expect any phase transitions when > 0. On 
the other hand, if is negative , so that <j> = 0 is a maximum at zero temperature, we 
can balance the negative zero-temperature mass against the positive contribution from 
the finite temperature piece. Thus, we expect that there will be a phase transition at 
some critical temperature Tc in this case. 

Whether a phase transition occurs depends upon the sign of m 2 as can be seen by 
examining A Vr(4>). An example of the temperature-dependent part of the potential is 
shown in Fig. 4. Clearly the curvature at (j> — 0 becomes more positive as the tem- 
perature increases. This does not depend upon the sign of the curvature of the zero- 
temperature potential as AVx(<£) is independent of mj. From Fig. 4 we also see that the 
finite-temperature corrections will always increase V(7r/2) more than V^(O), so if at zero 
temperature 7 r /2 is a maximum of the potential, it will remain so at high temperature. 
Of course, the actual value of m 2 is arbitrary, since it contains a renormalization coun- 
terterm. The value of ml is only technically naturally small, since it is protected by the 
chiral symmetry (and the residual discrete symmetries). 

There is no analytic expression for Tc] however, we can show that Tc must be of 
order m±. First we show that Tc cannot be much larger than the fermion masses m± 
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by means of the high temperature expansion. 27 For m/T <C 1 we can expand the finite 
temperature potential as: 

(5.12) 


AVrM = E i £j*J(*)ln^P + 


J=+,~ 




where we neglect terms such as (T-dependent) constants or (T-independent) parts de- 
pending on cos 2 <f>/ f. These terms are unimportant as far as computing the effective 
mass. The critical temperature is obtained by setting the second derivative of the full 
potential at <j> = 0 to zero. Doing this within the high temperature approximation yields 


2 mi-m 3 + ( 2 m\ 2 mi\ 


(5.13) 


We can solve this for Tc: 

~ rn \ ) 


=m+ fe) 


exp 


47 r 

X m 


ir 2 f 2 ml \ 
- m +) / 


(5.14) 


(recall < 0). It is easily seen that this quantity is at most of order m± so that the 
conditions for the validity of the high temperature expansion do not obtain, and the 
phase transition cannot occur at T m±. Now consider the possibility that Tc is much 
less than m±. In the limit mjT » 1, clearly A Vr(<f>) <x exp(— m±/T), so the phase 
transition cannot take place at T <C m±. It follows that we should expect the phase 
transition to occur near the scale set by m±. 

In Fig. 5 we show the total ^-dependent potential as a function of temperature. 
Clearly there is a phase transition somewhere in the range 3m_ < Tc < 5m_ when 
the high-temperature maxima become the low- temperature minima. Just as clearly, the 
phase transition will be second order. A unique feature of this model is that at the critical 
temperature the potential is absolutely flat — <f> becomes a free field (not simply massless 
as in a typical second-order transition). This can be understood by observing that the 
only extrema of the potential are at <f> = 0, 7t/2, it and when the full potential 
evaluated at 7t/2 becomes equivalent to the potential evaluated at 7r, there can be no 
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intervening extrema so the potential must be flat. Thus at the critical temperature the 
Z 2 symmetry is promoted to a (non-linearly realized) f7( 1) symmetry. 

Above the critical temperature there is still a Z 2 symmetry — there is in no sense 
a larger symmetry at high temperature. Another interesting feature of this model is 
that at high temperature, the potential becomes T independent (except for an additive 
T-dependent, ^-independent constant). This in in contrast to usual high-temperature 
scalar field theory where the mass of the scalar field at high temperature is proportioned 
to T. 

Notice that the phase transition can lead to the formation of domain walls. For 
instance, if <f> is at the minimum <f> = 7r for T > Tc, when the phase transition is complete 
regions of the Universe with <f> = 7r/2 will be separated from regions with <j> = 37r/2 by 
a domain wall. However there is one concern with the above scenario: There may be 
no physical mechanism to set <f> to its high-temperature minimum. The value of <f> at 
high temperature may be free to roam and may not be pinned to any particular value. 
This is because the Z 2 symmetry of the model implies that the <j> 2 T 2 term will not be 
present, and at high temperatures the leading temperature-dependent, ^-dependent term 
is ^ 3 ln(T 2 ), which grows slowly with T. One might well imagine a scenario where <f> has 
insufficient time to relax to its high-temperature minimum before the onset of the phase 
transition. If <j> has a value away from the low-temperature maximum at the onset of the 
phase transition, and it is constant throughout the Universe (say set during inflation), 
then the entire Universe may evolve to the same low-temperature value of <f> and domain 
walls would not appear. Therefore, if domain walls are produced, / must be considerably 
less than the scale of inflation. 

Finally, we digress for a moment to make sure we know just exactly whose temperature 
enters into the above expressions. Recall that the light neutrinos decouple from the 
ambient plasma at T& ~ 1 MeV. Thus after this time the neutrinos are not in thermal 
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equilibrium. However, the neutrino distribution function is still that of a particle in 
thermal equilibrium (so long as T is not too much less than the mass) with an effective 
temperature given by u(ti?)T£)/a(t) where &(£) is the scale factor, and tp is the cosmic 
time at which decoupling occurs. Thus, it is this effective temperature that appears in 
the finite-temperature effective potential. 

B. ^-symmetric models 


It is simple to generalize the models in the previous section, with its Z 2 symmetry 
amongst the fermions to one with N fermions and a corresponding Z n symmetry. The 
Yukawa couplings for such a model are: 

-£yuk = Y $3 { m + e[cos(^// + 2 *j/ N ) + *7s sin (<£// + 2irj/N)]}ipj. (5.15) 

3=0 

This theory has the Zn symmetry given by 

Mf->*/f + 2T/N, (5.16) 


where now the index i is taken mod N. 

The same methods used in the Zj case can be used here to calculate the effective 
potential. We find [cf., Eq. 5.3)] 

v.w = | [£ - 1 )\ • < 517 > 

with A being the ultraviolet cutoff, as usual, and 


M?((f>) = m 2 -f e 2 + 2mecos(<£// + 2 irj/N) 


(5.18) 


for j = 0, • • • , N — 1. In parallel with the Zj case we introduce an arbitrary scale /i, and 
rewrite the potential as [cf., Eq. (5.4)] 

N ~ i r a 4 a 2 

VJS) = Y — 

^ 32tt 2 8tt 2 lKY> 
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(5.19) 


1 

167T 2 






+ 


1 

16tt 2 



Our next task is to ascertain what types of counterterms must be included to absorb 
the divergences present in (</>). Before proceeding, we note the following: 

= N{m 2 + e 2 ) (for fV > 1) 

i= o 

= N[(m 2 + e 3 ) 2 + 2m J e 3 ] (for JV > 2). (5.20) 

i=o 

Whereas in the Z 2 case £j_ 0 M 2 ((f>) was <f> independent, if the discrete symmetry is Z^> 2 , 
then X!]=o M*(4>) is also <j> independent. Thus, the only counterterm we need to add is 
the ^-independent term Vo, given by 


v. - VM - E (£ - £*?<« + ' 


(5.21) 


32ir a 8x r 

For ming the total potential and dropping irrelevant ^-independent terms, we find 

(5.22) 


N-l 


VM = VnM - £ 


Again V 0 can be found by some renormalization condition, and the fi dependence in Vo(/z) 
will cancel the /z dependence in the log term rendering the entire potential finite and /z 
independent. Since we did not need to add any ^-dependent counterterms, the <j> mass 
is calculable in terms of the parameters of the theory, i.e., m and e. 

The extrema of are somewhat trickier to find than in the Z 2 case. It can, 

however, be shown that these are located at <)>// = 0, ir/N mod 2irjN. Whether these 
are maxima or minima depends on N\ for N even (odd) <f> = 0 is a max (min) while 
(f>j f = tt/N is a min (max). The potential again has a simple periodic form. The form 
of the potential for JV = 3 is shown in Fig. 6. 

Now the temperature corrections are easy to calculate — they are given by Eq. (5.10) 
where now the sum on j runs from 0 to N — 1. An example of the temperature-dependent 
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corrections to the potential is shown in Fig. 7. It is clear that the sign of the temperature- 
dependent part of the potential is opposite to the sign of the zero-temperature potential. 

The total potential V{<f>) + AV T (<j>) is shown in Fig. 8 for several temperatures. The 
interesting result is that at high temperatures the Z N symmetry is promoted to an 
exact (non-linearly realized) 17(1) symmetry. That this should occur is easy to see by 
examination of the high-temperature expansion of the finite-temperature potential [cf., 
Eq. (5.12)]: 

N ~ 1 1 M?(6) 

AW)= (5.23) 

which exactly cancels the entire ^-dependent part of Eq. (5.22). Thus at high temperature 
the potential becomes exactly flat. 

What cosmology might one expect given the temperature behavior of the potential? 
Clearly between T ~ / (when the effective potential makes sense) and T ~ m, <j> is free 
to take on any value. Below some temperature of order m the potential minima will start 
to become important and different regions of the Universe will have different values of <f> 
with domain walls between them. Thus, effectively there is a phase transition at T ~ m 
where the order parameter (in this case (f > ) evolves from whatever value it had at high 
temperatures to a zero-temperature minimum. In this case the transition is similar to 
the phase transition associated with axions, although we emphasize that the underlying 
dynamics Me quite different in the two cases. 

Also in analogy with the axion case, if inflation occurs at a scale less than /, then 
one might expect (f> to be set to a single value throughout the Universe. If this happens, 
when the transition occurs there will be a single initial value of <j> that will be random, 
there is nothing to perch the initial value of <f> on a low-temperature maximum, and the 
Universe will most likely end up in a single value of 4 > — no domain walls. 
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VI. CONCLUSIONS 


In this paper we have given a general discussion of the thermal physics of pseudo- 
Nambu-Goldstone bosons. These afford a natural way of generating very soft scales of 
potential interest to astrophysics, typically of order ~ where we might 

choose m, ma n ~ m neut rino ~ 10 -2 eV, and /urge ~ f gut ~ 10 15 GeV. These objects have 
a precedent in elementary particle physics in the familiar nucleon-meson system, as well 
as scores of theoretical generalizations, and the models considered here really involve 
no additional physical components. As in the case of the invisible axion, or familons, 28 
arions, etc., we are simply abstracting the scales to those that are of potential interest 
to astrophysics or cosmology. We have discovered that the thermal behavior of these 
systems is very simple, controlled largely by the residual symmetries of the low energy 
potential. 

Though we have largely focused on the specific models of Ref. (7), these models 
capture most of the physics that can generally occur in the context of pNGB’s. For 
example, the Zff models for large N have very soft breaking of the continuous 1/(1) 
symmetry, due to the cancellation of the fermion loops at high momentum from the 
discrete symmetry. It is, therefore, not surprising that the thermal behavior of this 
system imitates that of the axion, since the PQ symmetry of the axion is broken only in 
the far infrared limit of QCD. 

To a good approximation we may summarize the thermal physics as follows. The 
^-dependent part of the potential has the form: 

V{<j > ) = c(T)m 4 cos(N<j>/f) (6*1) 

where c(T) is a smoothly varying function of T with the following possible behaviors: 

1. c(T) is slowly varying with no sign change over the full range of temperatures 
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(o < T < /); 


2. c(T ) is slowly varying with a sign change for T ~ m, as in the (Z 2 model); 

3. c(T) is slowly varying with asymptotic zero c(T) — + 0 as T — ► oo (as in Zpr, for 
N > 2). 

Here we have not addressed the issue of cosmological implications (if any) for the 
formation of structure or other possibile signatures. It seems that the options here are (i) 
to pursue schemes that lead to soft domain walls, or other topological configurations, with 
“thicknesses” of order m^ 1 , which form after the 3°K microwave background decoupling, 8 
or (ii) to try to buld a natural version of the Press, Ryden, Spergel 11 scheme. 29 The latter 
has an additional potential fine-tuning problem associated with initial conditions that 
may be remedied in something like the Z 2 scheme with a sign change in c(T). We also 
mention that other large scale signatures, such as periodic redshifts, might require some 
bizarre version of schemes as discussed here. 30 

If the symmetries and dynamics of particle physics are a guide, then it seems likely 
that either ultra-low mass fermions, such as massive neutrinos, or ultra-low mass bosons, 
such as pNGB’s, are the best candidates for potential new cosmological effects. Re- 
stricting attention to such classes of particles is a powerful simplification rather than a 
complication. The existence of such objects implies dramatic new physics at the high- 
est energies, O(f), that lead to phenomena on the largest distance scales (as large as 
//m 2 ), which are of relevance for cosmology. We thus feel that the general discussion 
of the thermal behavior of pNGB’s given here is an important consideration for future 
cosmological model building efforts. Cosmologists should learn the physics of pNGB’s 
and think about their potential implications in the early and not-so-early Universe. 
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FIGURE CAPTIONS 


Figure 1: Tadpole diagrams used in calculating the effective potential. Dashed lines rep- 
resent bosons and solid lines represent fermions. 

Figure 2: Loops that lead to induced pNGB mass terms in (a) a chiral Lagrangian scheme 
of Eq. (4.8) and (b) the Majorana scheme of Eq. (4.15). 

Figure 3: The zero-temperature potential of the Z 2 model. 

Figure 4: The temperature- dependent corrections to the Z 2 model. 

Figure 5: The total potential of the Z 2 model. 

Figure 6: The zero-temperature potential of the Z$ model. 

Figure 7: The temperature-dependent corrections to the Z$ model. 

Figure 8: The total potential of the Z 2 model. 
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